function [Prob_H_M_L,c]=LRmodel(termN,ddG3)


%parameters:
AH=0.2973;AM=0.1661;
BH(1)=1.8767; BM(1)=0.9852;
BH(2)=0.3737; BM(2)=0.0193;
BH(3)=1.8744; BM(3)=1.3088;
BH(4)=1.7387; BM(4)=1.0301;
BH(5)=0.4764; BM(5)=0.0406;
BH(6)=1.5084; BM(6)=1.0877;
BH(7)=-0.2463; BM(7)=-0.4025;
BH(8)=-0.6404; BM(8)=-0.3337;
BH(9)=0.2381; BM(9)=-0.081;
BH(10)=-0.7873; BM(10)=-0.2945;
BH(11)=-0.5707; BM(11)=-0.3471;
BH(12)=-0.3552; BM(12)=0.0987;
BH(13)=-1.1939; BM(13)=-0.6077;
BH(14)=-2.2329; BM(14)=-1.3262;
BH(15)=-1.1755; BM(15)=-0.9055;
BH(16)=-2.2846; BM(16)=-1.2957;
CH=-0.0498; CM=0.2145;

N=size(termN,1);
Prob_H_M_L=zeros(N,3);
c=zeros(N,1);
for i=1:N
    highE=exp(CH+AH*ddG3(i,1)+BH(termN(i,1)));
    medE=exp(CM+AM*ddG3(i,1)+BM(termN(i,1)));
    totalE=1+highE+medE;
    Prob_H_M_L(i,1)=highE/totalE;
    Prob_H_M_L(i,2)=medE/totalE;
    Prob_H_M_L(i,3)=1-Prob_H_M_L(i,1)-Prob_H_M_L(i,2);
    [temp,c(i)]=max(Prob_H_M_L(i,:));
    
end
return;